A RECOMBINATION ALGORITHM FOR THE 
DECOMPOSITION OF MULTIVARIATE RATIONAL FUNCTIONS 
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Abstract. In this paper we show how we can compute in a deterministic way 
the decomposition of a multivariate rational function with a recombination 
strategy. The key point of our recombination strategy is the used of Darboux 
polynomials. We study the complexity of this strategy and we show that this 
method improves the previous ones. In appendix, we explain how the strategy 
proposed recently by J. Berthomieu and G. Lecerf for the sparse factorization 
can be used in the decomposition setting. Then we deduce a decomposition 
algorithm in the sparse bivariate case and we give its complexity. 



Introduction 

The decomposition of univariate polynomials has been widely studied since 1922, 
see [Rit22| . and efficient algorithms are known, see |AT85| IBZ85| IKL89| IGat90a| 
IGat90b[ IGie88[ IKlu99] . There also exist results and algorithms in the multivariate 

case |Dic87l TGat90a[ IGie88l IGGR03] . 

The decomposition of rational functions has also been studied, |Gie88} |Zip91| 
IAGR95[ IGW95] . In the multivariate case the situation is the following: 
Let f{X 1 ,...,X n ) = f 1 (X 1 ,...,X n )/f 2 (X u ...,X n ) £ K(X u ...,X n ) be a ratio- 
nal function, where K is a field and n > 2. It is commonly said to be composite if it 
can be written / = u(h) where h(Xi, . . . , X n ) £ K(Xl, . . . , X n ) and u £ K(T) such 
that deg(u) > 2 (recall that the degree of a rational function is the maximum of 
the degrees of its numerator and denominator after reduction), otherwise / is said 
to be non-composite. In this paper, we give an algorithm which computes a non- 
composite rational function h £ K(Xi, . . . , X n ) and a rational function u £ K(T) 
such that / = u(h). 

In [ChelOj , the author shows that we can reduce the decomposition problem to a 
factorization problem and gives a probabilistic and a deterministic algorithm. The 
probabilistic algorithm is nearly optimal: it performs 0(d n ) arithmetic operations. 
The deterministic one computes 0(d 2 ) absolute factorizations and then performs 
0(d n+u+2 ) arithmetic operations, where d is the degree of / and ui is the feasible 
matrix multiplication exponent as defined in [GG03, Chapter 12]. We recall that 
2 < lj < 2.376. As in ChelS], we suppose in this work that d tends to infinity and 
n is fixed. We use the classical O and O ("soft 0") notation in the neighborhood 
of infinity as defined in |GG03[ Chapter 25.7]. Informally speaking, "soft 0"s are 
used for readability in order to hide logarithmic factors in complexity estimates. 
In this paper we improve the complexity of the deterministic algorithm. With 
this algorithm we just compute two factorizations in K[Xi, . . . ,X n ] and then we 
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use a recombination strategy. Under some hypotheses this new method performs 
(D(d n+U1 ~ 1 ) arithmetic operations. 

The decomposition of multivariate rational functions appears when we study the 
kernel of a derivation, see }MO04j . In |MQ04| the author uses Darboux polynomials 
and gives an algorithm which works with 0{d un ) arithmetic operations. 
In this paper, we are also going to use Darboux polynomials (see Section Q] for a 
definition) and we add a recombination strategy. Roughly speaking, we are going 
to factorize the numerator and the denominator and then thanks to a property of 
Darboux polynomials we are going to show that we can recombine the factors and 
deduce the decomposition. 

The decomposition of multivariate rational functions also appears when we study 
intermediate fields of an unirational field and the extended Liiroth's Theorem, see 
GRSOlj lChelO] . and when we study the spectrum of a rational function, see [Che 10] 
and the references therein. 

The study of decomposition is an active area of research: for a study on mul- 
tivariate polynomial systems see e.g. [FP091 IFGPlOj . for a study on symbolic 
polynomials see e.g. |Wat09j . for a study on Laurent polynomials see e.g. |Wat08] . 
for effective results on the reduction modulo a prime number of a non-composite 
polynomial or a rational function see e.g. |CN101 [BDN09, BCN], for combinatorial 
results see e.g. |Gat08j . 

In this paper, we improve the strategy proposed in MO04]. As in |MO04) . we 
consider fields with characteristic zero. Furthermore, as we want to give precise 
complexity estimate we are going to suppose that: 

Hypothesis (C): 

K is a number field: K = Q[a], a is an algebraic number of degree r. 
As in IChelOj . we are going to suppose that the following hypothesis is satisfied: 

Hypothesis (H): 

{(i) deg(/i -I- A/2) = deg Xn (fi + A/ 2 ), where A is a new variable, 
(«) R(A) = Res Xn (/x(0, X n ) + A/ 2 (0, X n ), d x Ji(0, X n ) + AB Xn M X n j) + 

where deg^-^ / represents the partial degree of / in the variable X n , deg/ is the 
total degree of / and Resx n denotes the resultant relatively to the variable X n . 

This hypothesis is necessary, because we will use the factorization algorithms 
proposed in [Lec07] . where this kind of hypothesis is needed. Actually, in |Lec07] 
the author studies the factorization of a polynomial F and uses hypothesis (L), 
where (L) is the following: 

Hypothesis (L): 

{(i) deg Xn F = deg F, and F is monic in X n , 
(ii)Res Xn (F(0,X n ),-^{0,X n )) ^ 0. 

If F is squarefree, then hypothesis (L) is not restrictive since it can be assured by 
means of a generic linear change of variables, but we will not discuss this question 
here (for a complete treatment in the bivariate case, see |CLQ71 Proposition 1]). 
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Roughly speaking, our hypothesis (H) is the hypothesis (L) applied to the poly- 
nomial fi + A/2. In (H,z) we do not assume that f\ + A/2 is monic in X n . Indeed, 
the leading coefficient relatively to X n can be written: a + Ab, with a,J £ K. In 
our algorithm, we evaluate A to A such that a + Xb ^ 0. Then we can consider 
the monic part of f\ + A/2 and we get a polynomial satisfying (L,i). Then (H,i) 
is sufficient in our situation. Furthermore, in this paper, we assume /1//2 to be 
reduced, i.e. f\ and fi are coprime. We recall in Lemma [9] that in this situation 
/1 + A/2 is squarefree. Thus hypothesis (H) is not restrictive. 

Furthermore, hypothesis (H) will also be useful in a preprocessing step, see Sec- 
tion[2] In this preprocessing step we reduce the decomposition to two factorizations 
of squarefree polynomials. 

Complexity model. In this paper the complexity estimates charge a constant 
cost for each arithmetic operation (+, — , x, -=-) and the equality test. All the 
constants in the base fields are thought to be freely at our disposal. 

In this paper we suppose that the number of variables n is fixed and that the 
degree d tends to infinity. 

Polynomials are represented by dense vectors of their coefficients in the usual 
monomial basis. For each integer d, we assume that we are given a computation 
tree that computes the product of two univariate polynomials of degree at most 
d with at most 0(d) operations, independently of the base ring, see GG03, The- 
orem 8.23]. Then with a Kronecker substitution we can compute the product of 
two multivariate polynomials with degree d with n variables with 0{d n ) arithmetic 
operations. We also recall, see |GG03[ Corollary 11.8], that if IK is an algebraic 
extension of Q of degree r then each field operation in K takes 0(r) arithmetic 
operations in Q. 

We use the constant lu to denote a feasible matrix multiplication exponent as defined 
in }GG 03, Chapter 12]: two n x n matrices over K can be multiplied with 0{n u ) 
field operations. As in [BP94j we require that 2 < ui < 2.376. We recall that the 
computation of a solution basis of a linear system with m equations and d < m un- 
knowns over K takes e^mcT -1 ) operations in K [BP941 Chapter 2] (see also [StoOO, 
Theorem 2.10]). 

In |Lec06| ILec07] the author gives a deterministic algorithm for the multivariate ra- 
tional factorization. The rational factorization of a polynomial / is the factorization 
in K[X], where K is the coefficient field of /. This algorithm uses one factorization 
of a univariate polynomial of degree d and CD(d n+u ^ 1 ) arithmetic operations, where 
d is the total degree of the polynomial and n > 2 is the number of variables. 

Main Theorem. The following theorem gives the complexity result of our algo- 
rithm. 

Theorem 1. Let f be a multivariate rational function in Q[a](Xi, . . . , X n ) of 
degree d, where a is an algebraic number of degree r. Under hypotheses (C) and (H), 
we can compute in a deterministic way the decomposition of f with (5(r(i n+w_1 ) 
arithmetic operations over Q plus two factorizations of univariate polynomials of 
degree d with coefficients in Q [a] . 

Comparison with other algorithms. There already exist several algorithms for 
the decomposition of rational functions. They all use the same global strategy: first 
compute h, and then deduce u. The first step is the difficult part of the problem. 
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In [Che 10] . we explain how we can perform the second step, i.e. compute u from h 
and /, with 0(d n ) arithmetics operations. 

In [GRSQlj , the authors provide two algorithms to decompose a multivariate ratio- 
nal function. These algorithms run in exponential time in the worst case. In the first 
one we have to factorize polynomials with 2n variables fi (X) f?.(Y) — fi (Y) f?.(X) 
and to look for factors of the following kind /ti (X)h?(Y) — hi (Y)h?(X). The au- 
thors say that in the worst case the number of candidates to be tested is exponential 
in d = deg(/i//2). Indeed, the authors test all the possible factors. 
In the second algorithm, for each pair of factors {h\,h 2 ) of /1 and f% (i.e. h\ di- 
vides /1 and h 2 divides /2), we have to test if there exists u £ K(T) such that 
/1//2 = u(h\/h,2). Thus in the worst case we also have an exponential number of 
candidates to be tested. 

To the author's knowledge, the first polynomial time algorithm is due to J. Moulin- 
Ollagnier, see |MO04] . This algorithm relies on the study of the kernel of the 
following derivation: S U1 (F) — u) A dF, where F G K[X] and lo = f 2 dfi — fidf 2 - In 
MO04] the author shows that we can reduce the decomposition of a rational func- 
tion to linear algebra. The bottleneck of this algorithm is the computation of the 
kernel of a matrix. The size of this matrix is 0(d n ) x 0(d n ), then the complexity 
of this deterministic algorithm belongs to 0(d nu ). 

The reduction of the decomposition problem to a factorization problem is clas- 
sical, see e.g. |Klii99l IGie88l IGW951 IGRSOlj . In jGhelOj the author shows that 
if we choose a probabilistic approach then two factorizations in K[Xi, . . . , X„] arc 
sufficient to get h and furthermore we do not have a recombination problem. This 
gives a nearly optimal algorithm. For the deterministic approach the author uses a 
property on the pencil /1 — A/2 and shows that with 0(d 2 ) absolute factorization 
(i.e. factorization in the algebraic closure of K) we can get h. This deterministic 
strategy works with 0(d n+u+2 ) arithmetic operations. 

In this paper, we are going to show that we can obtain a deterministic algorithm 
with just two factorizations in K[Xi, . . . ,X n ] and a recombination strategy. Our 
algorithm uses at most 0(d n+u ~ 1 ) arithmetic operations. This cost corresponds to 
the cost of the factorization and the recombination step. 

Our recombination problem comes from this factorization: 
If /1//2 = ui/u 2 (hi/h 2 ) then 

/1 - A/ 2 = e(hi - tih 2 ) ■■■(hi- t k h 2 ) 

where A, e 6 K, k = deg(ui/u 2 ) and t\ are the roots of the univariate polynomial 
mi(T) - \u 2 (T), see Lemma [TO] 

Thus with the factors hi — t\h 2 and h\ — t 2 h 2 we can deduce h. Unfortunately these 
factors are not necessarily in K[Xi, . . . , X n ] and are not necessarily irreducible. In 
this paper we show how we can reduce the problem to a factorization problem in 
K[Xi, . . . , X n ] and how we can recombine the irreducible factors of /1 — A/2 to get h. 

We can see our recombination scheme as a logarithmic derivative method. Roughly 
speaking, the logarithmic derivative method works as follow: 

If F(X, Y) = n* = i Fj(X, Y), where Jj-(X, Y) E A and A D K[X, Y] (for example 
A = K[[X]][V]), then we can write the irreducible factors F^X, Y) £ K[X, Y] of F 
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in the following way: Fi — Yij=i ^ ' 3 > where e,,j G {0, 1}. Thus we just have to 
compute the exponents e^j to deduce Fj. We compute these exponents thanks to 



this relation: 



3 = 1 ^ J 

With this relation the exponents e, j are now coefficients, and we can compute them 
with linear algebra. 

This strategy has already been used by several authors in order to factorize poly- 
nomials see e.g. [BHKS091 IBLS+041 ILecOSI RTCD71 rWeiTDj . Here, we use this kind 
of technique for the decomposition problem. With this strategy the recombination 
part of our algorithm corresponds to the computation of the kernel of a 0(d n ) x 0(d) 
matrix. 

In our context, we do not use exactly a logarithmic derivative. We use a more 
general derivation, but we use the same idea: if a mathematical object transforms 
a product into a sum then the recombination problem becomes a linear algebra 
problem. In this paper this mathematical object is the cofactor, see Proposition |8] 

Structure of this paper. In Section 1, we recall some results about the Jacobian 
derivative and Darboux polynomials. In Section 2, we describe a reduction step 
which eases the recombination strategy. In other words we explain how we can 
reduce the decomposition problem to a factorization problem. In Section 3, we 
explain how we can get h with a recombination strategy. In Section |4j we describe 
our algorithm with two examples. In Section [5] we conclude this paper with a 
remark on Darboux method and the logarithmic derivative method. In appendix, 
we explain how the strategy proposed recently by J. Berthomieu and G. Lecerf in 
BL10 for the sparse factorization can be used in the decomposition setting. Then 
we deduce a decomposition algorithm in the sparse bivariate case and we give its 
complexity. 

Notations. All the rational functions are supposed to be reduced. 
Given a polynomial /, deg(/) denotes its total degree. 

Given a rational function / = /1//2, deg(/) denotes max (deg(/i), deg(/2)). 

For the sake of simplicity, sometimes we write K\X] instead of . . . , X n ], for 

n > 2. 

u o h means u{h). 

Res(A, B) denotes the resultant of two univariate polynomials A and B. 
\S\ is the cardinal of the set S. 



1. Derivation and Darboux polynomials 

We introduce the main tool of our algorithm. 

Definition 2. A K-derivation D of the polynomial ring K[Xi, . . . , X n ] is a K-linear 
map from . . . , X n ] to itself that satisfies the Leibniz rule for the product 

D(f.g)=D(f).g + f.D(g). 

A K-derivation has a unique extension to K(Ai, . . . , X n ) and then we will also 
denote by D the extended derivation. 
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Definition 3. Given a rational function /1//2, the Jacobian derivative associated 
to /1//2 is the following vector derivation, i.e. an (n — l)-tuple of derivations: 



n-1 



D fl/f2 :K[X 1 ,...,X n ] — > (K[X 1 ,...,X n ]j 

(d Xl (hlh)dx 2 F - d X2 {h lh)d Xl 



F — > fl 

\d Xl {hlh)d Xn F - d Xn (fi/f 2 )d Xl Fj 

The Jacobian derivative has the following property: 

Proposition 4. Given f = /1//2 and g £ K(Xi, . . . , X n ) \ K the following propo- 
sitions are equivalent: 

(1) The rank of the Jacobian matrix 



Jac{f,g) 



I 'If df_ 

dXi dX, 

dg dg 

\dX 1 "' dX, 



is equal to one; 

(2) %//,(<?) = 0; 

(3) there exists h in K.(X±, , . . , X n ) such that f — u(h) and g — v(h) for 
u,v e K(T). 

Proof. See |PI07j for a proof. In |PI07j . K is supposed to be algebraically closed. 
However, we can remove this hypothesis because we have the equivalence: / is 
composite over K if and only if / is composite over K, see e.g. [BCN[ Theorem 
13]. □ 

Definition 5. Given D a vector derivation i.e. an m-tuple of derivations, a polyno- 
mial F € K[X] is said to be a Darboux polynomial of D if there exists Q € (K\X]j 
such that D(F) — F.Q. Q is called the cofactor of F for the derivation D. 

We deduce easily the following classical propositions. 

Proposition 6. f\ and fy. are Darboux polynomials ofDf ± /f 2 . 

Proposition 7. D t x i / 2 (/ii//i2) = if and only if hi and hi are Darboux polyno- 
mials with the same cofactor. 

The following proposition is the main tool of our algorithm. Indeed, this propo- 
sition shows that cofactors transform a product into a sum. Then thanks to the 
cofactors it will be possible to apply a kind of logarithmic derivative recombination 
scheme. 

Proposition 8. Let F £ K[Xt, . . . , X n ] be a polynomial and let F — F^ 1 ■ ■ ■ Ffj 

be its irreducible factorization in K[X 1; . . . , X n ]. Then: 

F is a Darboux polynomial with cofactor Gf if and only if all the Fi are Darboux 
polynomials with cofactor Qfi- Furthermore. Qf = &\Qf x + ' ' ' + &rQF r - 

Proof. See for example Lemma 8.3 page 216 in [PL A06j . □ 
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2. Reduction to a rational factorization problem 

In this section, we recall how the decomposition problem can be reduced to a 
factorization problem. Furthermore, we show that we can reduce our problem to a 
situation where f\ and f 2 are squarefree. First, we recall some useful lemmas. 

Lemma 9. // /1//2 is reduced in K(Xi, . . . , X n ), where n > 1 and A is a variable, 
then fi + A/2 is squarefree. 

Lemma 10. Let h = h\jh<x be a rational function in 'K(X). u — u\/u 2 a rational 
function in K(T) and set f = u o h with f = /1//2 G 'K(X). For all A G K such 
that deg(ui — AW2) = degu, we have 

fi - A/ 2 = e(/i! - t 1 h 2 ) ■■■(hi- t k h 2 ) 

where e G K, k — degu and ti are the roots of the univariate polynomial u\(T) — 
Xu 2 {T). 

Proof. See |ChelO| Lemma 8, Lemma 39]. □ 

Remark 11. If A = fi(a)/f 2 (a), where a = (ai, . . . , a n ) G K™, then we can suppose 
that ti G K. Indeed, ti = hi (a)/h 2 (a) G K. 

The following lemma says that we can always suppose that degui = deg«2 = 
degu. 

Lemma 12. Let h = h\/h 2 be a rational function in K.(X), u — u\/u 2 a rational 
function in K(T) and set f — u o h with f = /1//2 G K(X). There exists an 
homography H(T) = (oT + b)/(aT + 13) G K(T) suc/i tfiai; 

u o H = u\/u 2 , deg ui = deg M2, awe? / = — o h, where h = H^ 1 o ft, and _ff _1 is 

u 2 

the inverse of H for the composition. 

Proof. If deg Mi = degU2 then we set H(T) = T. 
If degit2 > deg ui then we have: 

H(H(T)) = IlS" 1 (aT + b-XjjaT + P)) , )degtt2 _ degUl 
U2 \ >) H^(aT + b- fk (aT + f3))' y 

where Wi(Aj) = and u 2 (/j,i) = 0. 
We set: 

dcg ui 

«i(T) - (aT + /3) dcgu ^ dcgtl1 . J| (aT + 6- A^aT + Z?)) 

i=l 

= ui(iJ(T)).(aT + /3) dog " 2 G K[T] 

degu 2 

ua(T) = f[ (aT + b- Ln(aT + p)) 

i=l 

= u 2 ( J ff(T)).(ar + /3) dcsU2 G K[T]. 

If a — A;a / 0, a ^ 0, and a — then we get degui = degU2 = degM2- 

To conclude the proof we just have to remark that degH = 1, thus H is invertible 

for the composition. □ 
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In order to ease the recombination scheme we reduce our problem to a situation 
where the rational function is squarefree, i.e. the numerator and the denominator 
are squarefree. The following algorithm shows that if f\ or f 2 are not squarefree 
then we can compute an homography U(T) £ K(T) such that U(fi/f 2 ) is square- 
free. Furthermore, if we know a decomposition U(f\/ f 2 ) = u(h) then we can easily 
deduce a decomposition /1//2 = U^ 1 (u(h)) . We recall that U is invertible for the 
composition because degU = 1. Now, we describe an algorithm which computes a 
good homography. 

Good homography 

Input: / = /1//2 6 K(Xi, . . . , X n ) of degree d, such that (C) and (H) are satisfied 
and a finite subset S of K n such that \S\ = 2d 2 + 2d. 

Output: U(T) = (T-X a )/(T-X b ) such that U(f) is squarefree, A a = fi/f 2 (a),X b = 
/1//2® where a,b e K n , X a ^ A 6 , and deg Xn (f 1 - X a f 2 ) = deg^J/i - A fc / 2 ) = d. 

(1) Compute J 1 (X n ) := fi(0,X n ), and J 2 (X n ) := / 2 (0, X n ). 

(2) Construct an empty list L. 

(3) For i from 1 to 2d 2 + 2d do: 

(a) Compute / := J 1 (i)/J 2 (i)y 

(b) If / ^ L then £ := concatenate^, [/]). 

(4) Construct an empty list C. 

(5) For k from 1 to 2d + 2 do: 

(a) Computed := Res Xn {j^X^-Llk]]^), d x J 1 {X n )-L[k]dxJ 2 {X n )). 

(b) If R '.^ and deg x (f 1 — L[k]f 2 ) = d, then C := concatenate(£, 

(6) A a :=C[l], X b :=£[2}. n 

(7) Return U{T) = (T - X a )/(T - X b ). 

Proposition 13. The algorithm Good homography is correct. 

Proof. In Step [3] we construct a list with at least 2d + 2 distinct elements because 
deg(/)=d. 

By hypothesis (H), i?(A) ^ and by |CG03I Theorem 6.22], deg(i?) < 2d- 1. Thus 
£ contains at least two distincts elements. 

As R(X a ) and R(Xb) are not equal to zero, and thanks to Step I5bl the condition on 
the degree is satisfied, we deduce that f\ — X a f 2 and f\ — Xbf 2 are squarefree. □ 

Proposition 14. The algorithm Good homography can be performed with at most 
0(d n ) arithmetic operations over K. 

Proof. Step[T]can be done with 0(d n ) arithmetic operations with Horner's method. 
In Step we use a fast multipoint evaluation strategy, then we can perform this 
step with at most 0(d 2 ) arithmetic operations, see |GG03[ Corollary 10.8]. 
In Step [5l the computation of the resultant can be done with €>{d) arithmetic 
operations, see [GG03, Corollary 11.16]. Thus Step [5] can be done with (D(d 2 ) 
arithmetic operations. 

In conclusion the algorithm can be performed with the desired complexity. □ 

Remark 15. Suppose /1//2 = vi/v 2 (h). With the algorithm Good homography we 
can write U(fi/f 2 ) — ui/u 2 {h) with u\ju 2 £ K(T), h e K(X\, . . . , X n ), and u\ 
(resp. u 2 ) has a root a\ (resp. a 2 ) in K. Indeed, we have u\ = v\ — X a v 2 (resp. 
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u 2 = vi — \bV2) and A a = /i//2(a) (resp. A;, = /1//2©) then we deduce that 
Q'i = hi/h 2 (a) (resp. a 2 = hi/h 2 (b)). 

3. The recombination method 

In this section we describe our recombination method. First, we introduce some 
notations. By Proposition HI F\ and F 2 are Darboux polynomials of Dp 1 ^p 2 .We 
denote by 

the cofactor of F k , where k = 1, 2, and Qp € . . . , X n ]. We set: 

F k = ]jF k , j 

for k = 1,2, and 

g, fe , = (a^,...,^). 

In Q[a][Xi, . . . , X n ] polynomials are denoted in the following way: 

r-i 

r=J2 E a ^ a£ *i T1 GQ[o][Jfi,...,x n ], 

|r|<de=0 

where a is an algebraic number of degree r, r = (n, . . . , r„), |r| = n H + T n , 

and a e , T G Q. We set 

coef(V, ofXj) = a e , T . 
Now we define the linear system S: 

Sl S2 

S .= J2 ^ijcoef(gg . , a e 2T) - J a; 2> jCoef (gg . , a e X T ) = 0, 
i=i j=i 
where |r| < d, < e < r — 1, and 2 < I < n. 

We denote by ker 5 the kernel of this linear system, and we remark that 

Sl S2 

x = (xi,i,...,a;2, S2 ) £ ker 5 <^=> ^".r ,,0' f -^£2,^2,, = 0- 



3=1 j=i 



We define the following maps: 



TTi : K Sl+S2 — > K Sl 
(xi,i,...,x 2}S2 ) 1 — > (xi^,...,Xx >Sl ) 

tt 2 : K Sl+S2 — > K S2 

(x ll x,...,X 2 , S2 ) 1 ► (x 2 ,l,-.-,X2, S2 ) 

The following proposition will be the key of our algorithm: 

Proposition 16. Suppose that F\/F 2 6 K(Xi, . . . ,X n ) comes from the algorithm 
Good Homography and F\jF 2 — u(h) where h = h\jh 2 £ K(Xi, . . . ,X n ) is a non- 
composite reduced rational function and u — u\/u 2 G K(T) is a reduced rational 
function, with degiti = degU2- 

We denote by u k = u k,i the factorization of in K[T], where k = 1, 2. 

We denote by F k = YViLi Fk,j the factorization of F k in K[Xi, . , . , X n ], where 
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k = 1, 2. 

Then: 

(1) u kti {^-).hf SUk * fl ri s eK[X 1 ,...,X n ] ande k ,ij G {0,1}- 

Furthermore, if we set e kl i ■= (ejfe,i,i, ■ ■ • , *Aen 

</ie vectors e k j, * = 1 , . . . , ife , are orthogonal for the usal scalar product. 

(2) FFe nave e*,* £ 70s(ker<S). 

(3) {e fei i, . . . , e fcjtfc } is a basis ofn k (kerS). 

Proof. (1) By Lemma [TOl applied to F1/F2 (resp. F2/F1) with A = 0, we get 

Fk = u k (hi/h 2 ).hi e * Uk =f[u k , i (h 1 /h 2 ).hf BUh - i . 

i=i 

Then we deduce 

u k ,i{hi/h2).hi esUKi = H in K[Xx, . . .,X n ] 

3=1 

with e k ^j g {0, 1} because F k are squarefree. Furthermore, the vectors e ki i 
are orthogonal for the usual scalar product because F k are squarefree. 
(2) We show this item for k = 1, the case k = 2 can be proved in a similar way. 
As F1/F2 comes from the algorithm Good Homography and as explained in 
Remark 1151 we can suppose that: 

u k .i{T) = (T- afc), with a fc g K. 

The previous item allows us to write: 



Ml,! 



W (ir : /v; ) .(/*) ' (11;,/^ ) 

By Proposition 2] applied to (t~) > we S e ^ then: 

U 2,l 2 



/ F 1 /F 2 \ S2 e 2 _ 1J .dc S u 1 



Now, we recall that F k _j are Darboux polynomials, see Proposition [6] and 
Proposition [5] Then by Proposition [FJ we deduce 

si s 2 

y^e 1;i jg Fl- . - deg(tt 1;i ) e2,i,jgjr 3J =0. 
j=i j=i 

It follows (ei^i, . . . , ei i i iS1) deg(wi >i ).e 2 ,i,i, . . . , deg(ui,j).e 2 ,i,s 2 ) G kerS. Thus, 
ei,i G 7Ti(ker5). 

(3) The vectors e^i, . . . , e k ,t k are linearly independant because they are or- 
thogonal. We just have to prove that these vectors generate n k (kerS). 
Suppose that p — (pi, . . . , p Sl+S2 ) G ker5. First, we clear the denominators 
and we suppose that p g Z Sl+S2 instead of Q Sl+S2 . 

In a first time we explain the strategy of the proof for this item, and in a 
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second time we will detail the proof. 
We set 

77 rr Sl F Pj 



where T\,T<x G ~K[X] and F\/F 2 is a reduced rational function. 



Our goal is to get this kind of equality: 

./•-, ;/•;:, \\l .11^ (11: 1^ P 



* ii; 11; ,..(11: ;/;y 



where m„ fc G N, J = {(1, 1), . . . , {h, 1), (1, 2), . . . , (i 2) 2)}, / num C /, 

7den C I and J num H /den = 0- 

By the unicity of the factorization in irreducible factors we deduce: 
Klip) = ^ rn Uli ei ti ~ 2J m u lti e i,i) 

(1.1) S/num (i,l)£Iden 

^(p) = ^ m« 2 /2,i - ^ m„ 2 i e2,i. 

(1.2) G/„um (i,2)G/ d o„ 

We get: {e^i, . . . , ek,t k } generates 7Tfc(ker5), and this is the desired re- 
sult. 

Now we detail the proof with four steps: 
(a) We remark: 

£1 = ui (h) = nS"(^i zB-iM 

IT 1 „. V / i—rdcffU/ 



f 2 ' n-: s r^i-M2^ 2 

where pk,i are roots of Uk- 



(b) We have: 



.h 2 , with k G Z, m« G 



Indeed, as p G ker<S, we have 

si s 2 

Thus YVjLi F\ 3 j an( i Iljli ■^2j L+3 are D ar boux polynomials with the 
same cofactor. By Proposition [3 we deduce: 



. .7 = 1 f 2,j 



Then D Fl ip 2 {J : \/3 : 2) — and thus T\jT 2 = v\/v 2 (h) by Proposition 
|U We denote by Aj the roots of v\ and v 2 and we get the desired 
result. 
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(c) We claim: 

n nLi Ume^ U,(hi/h 2 )h d 2 c ^A muk - 

-^r = -, rm ! where m u , G N. 

72 nLin(, fc ) e , don («M(^i/^)^ CSMfc -) w 

Indeed, we have: for all j there exists 6(j) such that Xj = Hs(j)- 

(To prove this remark we suppose the converse: There exists jo such 

that Xj ^ (J-k,ij for k = 1, 2 and i = 1, . . . , deg u. 

By definition of J-^ and by step !3b[ there exists (ki, ji) such that -Ffeij'i 
and /ii — \j h,2 have a common factor in C[Jf], We call "P this common 
factor. 

By step !3al there exists {k 2 , 22) such that P is a factor of h\ — /i k2t i 2 h2- 
Thus /ii — Xj h<2 and ft-i — /ifc 2 ,i 2 /i2 have a common factor. As Aj ^ 
A t fc 2 ,i 2 we deduce that P divides hi and ft.2- This is absurd because 
h\/h2 is reduced.) 

Thus k = 0, and for all j there exists G {1,2} and such that 

u k(j)(^j) = 0. 

As v\,V2 G Kf? 1 ], by conjugation, we deduce that if A 3 and Xj> are 
roots of the same irreducible polynomial u k j G K[T] then = m_y' . 
We denote by m Ufc . this common value. 

This gives the claimed equality with I num PI /den = because Txj T2 
is reduced. 

(d) Now we can prove equality {£). 

T x IlLi n ( i,fe)G/ num {uk,i{hi/h 2 )h d 2 CSUK 



T2 nLin (l , fc ) e / don (^,(^/^)^ cstife ' 
ri/c=i ri(i,fc)e/ num (rij-i ^kj 3 ) 
rife=i ri(i,fc)eJd.n ( n,-=i ^ ^ i-j ) 

This gives the desired equality (£). 



-, by step|3"cl 



by the hrst item. 



□ 



Now we describe our recombination algorithm: 



Recombination for Decomposition 

Input: / = /1//2 G K(Xi, . . . , X n ), such that (C) and (H) are satisfied. 
Output: A decomposition of / if it exists, with / = uoh, u = Ui/u^ with deg it > 2, 
and h = h\/h,2 non-composite. 



(1) Compute F = F1/F2 := U(f) with the algorithm Good homography. 

(2) For k=l, 2, factorizc F k = F k ,i in K[X] with F ktl irreducible. 

(3) For each F k .i compute the corresponding cofactor QF kA ■— D Fl /p 2 (F k ^)/F k ^. 

(4) Build the system S and compute the basis in reduced row echelon form B± 
of 7Ti (ker S) and B2 of 112 (ker S) . 

(5) For k=l, 2, find v k = (v kA , . . . ,«*.,«„) G B k such that: 

YfiL\ v k,i degFfe,, = minxes,, Yf^Li w i degF fcii , where w := (wi, . . . , w Sk ). 

(6) For k=l, 2, compute H k := KT ■ 
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(7) Set H := Hx/H 2 . 

(8) Compute u such that u{H) = f, 

(9) Return H, and u. 

Proposition 17. The algorithm Recombination for Decomposition is correct. 

Proof. Consider F1/F2 :— U(f). As we want to decompose /1//2, we just have to 
decompose F1/F2, because degU — 1 and then U is invertible. 
As F1/F2 comes from the algorithm Good Homography we can suppose, see Remark 
I15[ that Uk,i(T) = (T — otk) with ctk G K, and k — 1,2. Furthermore, by Lemma 
[12] we can also suppose that degui = degU2- 

Then by Proposition [1^1 the basis Bk of ^(kexjS) are . . . , efe,t fc }- 

The vector e k ,i gives the polynomial Hk,i = YljLi F kT' 3 = u k,i(h)h 2 esUh,i . 
Furthermore degHk,i — Ylj'Li e k,i,j deg-Ffcj = deg u k ,% deg h. Thus in Step [5] 

min Wi deg F fcji = deg h, 
%—\ 

because this minimum is reached with e^i € Sfe. Hence Wfe in Step [6] gives = 
Uk,i(k)(h)h^ C6UkMk} with degw M ( fe ) = 1. 

It follows H = (h 1 -ah 2 )/(h 1 ~f3h 2 ) with a, j3 G K. Thus H = v(h) with deg v = 1, 
then the algorithm is correct. 

□ 

Proposition 18. The algorithm Recombination for Decomposition can be performed 
with 0(rd' i+ " _1 ) arithmetic operations over Q and two factorizations of univariate 
polynomials of degree d with coefficients in K. 

We recall that in our complexity analysis the number of variables is fixed and 
the degree d tends to infinity. 

Proof. Step Q] uses (D(d n ) arithmetic operations over K by Proposition [Til thus it 
uses (D(rd n ) arithmetic operations over Q. 

Step [2] uses 0(d n+u) ~ 1 ) arithmetic operations over K because we can use Lecerf's 
algorithm, see Lec07 . Thus we use 0{rd n+u '^ 1 ) arithmetic operations over Q and 
two factorizations of univariate polynomials of degree d with coefficients in K. 
In Step [31 we compute Dp 1 /p 2 (Fk t i), thus we perform 2(n — 1) multiplications of 
multivariate polynomials. We can do this with a fast multiplication technique, 
and then this computation costs 0{nrd n ) arithmetic operations over Q. Then we 
divide Dp l /p 2 [F^^) by Fk,i- We have to perform n — 1 exact divisions, thus with a 
Kronecker subsitution we reduce this problem to n — 1 univariate divisions, and the 
cost of one such division belongs then to (D(rd n ). As si and s 2 are smaller than d, 
Step [3J costs 0(nrd n+1 ) arithmetic operations over Q. 

Step[4]needs 0(nrd n d u ~ 1 ) arithmetic operations over Q with Storjohann's method, 
see [StoOO, Theorem 2.10]. Indeed, S has 0((n — l)rd n ) equations and s\ + s 2 
unknowns, thus at most 2d unknowns. 

Step [5] has a ncgligeable cost because dimQ -Kk (ker S) = tk is smaller than d and Sk 
is also smaller than d. 

In StepO we use a fast multiplication technique and we compute Hk with 0(rd n ) 
arithmetic operations over Q. 

Step M can be done with 0(rd n ) arithmetic operations over Q, see [ChelOj . 
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Thus the global cost of the algorithm belongs to 0(rd n+UJ 1 ) arithmetic operations 
over O. □ 



4. Examples 

In this section we show the behavior of the algorithm Recombination for Decom- 
position with two examples. We consider bivariate rational functions with rational 
coefficients. Thus hypothesis (C) is satisfied. 



4.1. / is non-composite. We set: 

/i = (l + X + Y 2 ) (X + Y) = X + X 2 + XY 2 + Y + YX + Y 3 , 

f 2 = A - (Y 2 -X - 1) (Y- 2X + 1) = -X 2 + 3XY 2 + 2Y + 2YX -Y 2 + 1 

Wc have deg(/i + A/ 2 ) = deg y (/i + A/ 2 ) = 3, and 

ifea y (/i(0,y)+A/ 2 (0,r), 9r/i(0,F)+A^/ 2 (0,F)) = -4-24A-92A 2 -64A 3 +8A 4 . 

Thus hypothesis (H) is satisfied. 

The algorithm Good homography gives: 

A Q = /i(0,0)// 2 (0,0) = and A b = /i(0, l)// 2 (0, 1) = 1. 

Then 



Fi 


= (l + X + Y 2 )(X^ 


■Y), 




= l + X + Y 2 , 




Fl,2 


= X + Y 




F 2 


= (Y 2 -X-l)(Y- 


21+1) 


F 2 ,i 


= Y 2 - X - 1 




F 2 ,2 


= Y-2X + 1 





The cofactors are: 

G Fl l = 3X 2 + 8YX 2 + 2X - 2YX + 7XY 2 - 1 + 3Y 2 - 6Y 3 -6Y 4 + 2Y 
G Fl2 = 3 X 2 + 8 YX 2 + AYX + 61 - 6Y 2 - AY + 3- 3 F 4 - 2Y 3 
Gf 21 = 3X 2 + 8YX 2 + XY 2 - 6YX + 2X - I - 2Y - 6Y 4 -11Y 2 -6Y 3 
Gf 22 = 3X 2 + 8YX 2 + 6XY 2 + 8YX + 6X -3Y 4 + 8Y 2 - 2Y 3 + 3 
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The linear system S is the following: 



- -1 


3 


-1 


3 " 


2 


6 


2 


6 


3 


3 


3 


3 


o 


o 


o 


o 


o 


o 


o 


o 


2 


-4 


-2 


o 


—2 


4 


—6 


8 


8 



8 



8 



8 



3 


-6 


-11 


8 


7 





1 


6 














-6 


-2 


-6 


-2 














. -6 


-3 


-6 


-3 . 



A basis of ker(5) is given by: {(—1, — 1, 1, 1)}. 
Then it follows that /1//2 is non-composite. 

4.2. / is composite. Here we set: 

h x = (l + X + Y 2 ) (X + Y) 
h 2 = h 1 -(Y 2 -X-l){Y-2X + l) 
ui = T.(T-l) 
u 2 = T 2 + l 
/1//2 = u 1 /u 2 (h 1 /h 2 ). 

We have constructed a composite rational function /1//2 and now we illustrate 
how our algorithm computes a decomposition. We can already remark that in the 
previous example we have shown that h\jh 2 is non-composite. 

In this situation the hypothesis (H) is satisfied and the algorithm Good Homog- 
raphy gives: 

A a = /i(0,0)// 2 (0,0) = and A 6 - /i(0, 2)// 2 (0,2) - 90/101. 
Then: 

F M = 1 + X + Y 2 
F ia = 2X — Y — l 
Fi, 3 = Y 2 -X-l 
F M = X + Y 

F 2 ,i = 2X 2 + 11X + 9 + 29XY + 29Y + 38XY 2 -9Y 2 + 11Y 3 
F 2 , 2 = HX 2 + X- 10- 19YX- 19Y -29XY 2 + 1QY 2 + Y 3 
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The basis in reduced row echelon form of 7ri(ker<S) (resp. 7T2(ker5)) is {(1, 0, 0, 1); (0, 1, 1, 0)} 
(resp. {(1,0); (0,1)}). 

Step[H]in the algorithm Recombination for Decomposition gives: v\ — (1,0,0, 1) and 

W2 = (l,0). 

Then we have H\ :— F± t i.Ft t 4 and H 2 := i*2,i- 

We remark that Hi = hi and that H 2 = ll/ii + 9h 2 . Then Hi/H 2 = w(hi/h 2 ), 
where w(T) — T/(11T+ 9). As hi/h 2 is non-composite and degw = 1, we get a 
correct output. 

5. Conclusion 

In conclusion, we summarize our algorithm with a "derivation point of view" . 
In order to decompose /1//2, we have computed with Darboux method a rational 
first integral of f/w/j with minimum degree. That is to say we have computed 
hi/h 2 e K(Xi, . . . ,X n ) such that Dj 1 /f 2 (hi/h 2 ) = and deg(hi/h 2 ) is minimum. 
In a general setting, Darboux method works as follows: If we want to compute a ra- 
tional first integral of a derivation D, first we compute all the Darboux polynomials 
Fi and their associated cofactors Qpt, second we solve the linear system 

i 

Then thanks to Proposition^ we deduce that J\ i FT* is a first integral, i.e. .D(]X FT* ) = 0. 

When we consider the derivation Df x / f 2 the computation of Darboux polynomi- 
als is reduced to the factorization of fi + A/2 . Thus this step can be done efficiently. 
In the general setting, we can also reduce the computation of Darboux polynomials 
to a factorization problem, sec [Chcll . 

During the second step, we compute the kernel of J^. eiGFi = 0. It is actually a 
recombination step. Indeed, this system explains how we have to recombine Fi in 
order to get a rational first integral. Furthermore, the cofactor Qp i = D(Fi)/Fi can 
be viewed as a logarithmic derivative. 

In conclusion, the recombination scheme used in this paper is called nowadays the 
logarithmic derivative method, but this method is Darboux original method. 

Appendix A. Convex-dense bivariate decomposition 

In this appendix we give complexity results for the decomposition of sparse bivari- 
ate rational functions. These results rely on a strategy proposed by J. Berthomicu 
and G. Lecerf in uBLlO] . 

Given a polynomial f(X, Y) E K[X, Y], its support is the set Sf of integer points 
(i;j) such that the monomial X l Y 3 appears in / with a non zero coefficient. The 
convex hull, in the real space K 2 of Sf is denoted by N(f) and called the Newton's 
polygon of /. We denote by \N(f)\ the number of integral points of N(f). We 
called \N(f)\ the convex-size of /. 

Roughly speaking, the transformation proposed in |BL10] consists in a monomial 
transformation that preserves the convex-size but decreases the dense size. The 
considered transformation T can be described in the following way: 

T = B o £, where 
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B(X i Y j ) = x i+bl Y 3+b2 , 6i,6 2 e z, 
£pry J ') = X aii+a2j Y a3l+a ^ , ai a 4 - a 2 a 3 = ±1. 
T can be defined on K[X, Y, X -1 , F _1 ], and we define: TQ^ j AjX^) = 

The transformation £ corresponds to the linear map: A where 



.4 



<2i a 2 

We denote by £ _1 the transformation corresponding to A -1 . 

If /(X,Y) e K[X,y], then £(/) e K^y,!" 1 ^- 1 ] and £(/) can be written 
£(/) = c L (f)X (f), where £ (/) G K[*,Y] andc £ (/) = X*yj e KLY, F, X" 1 , Y' 1 } 
Furthermore, we also have C(F 1 .F 2 ) — C(F 1 ).C(F 2 ). 

Let S be a finite subset of Z 2 . Set S is said to be normalized if it belongs to 
N 2 and if it contains at least one point in {0} x N, and also at least one point in 
N x {0}. For such a normalized set, we write d x (resp. d y ) for the largest abscissa 
(resp. ordinate) involved in S, so that the bounding rectangle is 1Z — [0, d x ] x [0, d y \. 
The following result is proved in [BL101 Theorem 2]: 

For any normalized finite subset S ofZ 2 , of cardinality a, convex-size n, and bound- 
ing rectangle [0,^] x [0,^], and dense size S = (d x + l)(d y + l), one can compute an 
affine map T = Bo £, with O(crlog 2 S) bit- operations, such that T(S) is normalized 
of dense size at most 9n. 

We are going to use this transformation in order to prove: 

Theorem 19. Let f 1 /f 2 (X,Y) e K(X,Y) such that deg(/i// 2 ) = d, 2V(/i) C N, 
N(f 2 ) C N and N is normalized. Then 

(1) IfWL is field with characteristic or at least d(d— 1) + 1 and (H) is satisfied, 
then there exists a probabilistic algorithm which computes the decomposition 
0//1//2 with at most 0{\N\ 1,5 ) operations in K and two factorizations of 
a univariate polynomial of degree at most 9|iV| over K. 

(2) // (C) and (H) are satisfied, then there exists a deterministic algorithm 
which computes the decomposition of /1//2 with at most 0(r.\N\( u+i " 2 ) 
operations over Q and two factorizations of an univariate polynomials of 
degree at most 9\N\ over Q[a]. 

Now, we explain how we use the transformation T in the decomposition setting. 

Proposition 20. If h/ f 2 = u(hi/h 2 ) then T(f 1 )/T(f 2 ) = u{C(h{)l C(h 2 )) . 
IfT{fi)/T{h) = u(H 1 /H 2 ) then fr/f 2 = u^- 1 ^)/^ 1 ^)) . 

Proof. We prove the first item, the second can be proved in a similar way. 

We have: ^- — J}^ 1 tlhlhlL^ where [i^ , are roots of Ufe. Then, 

H HjC^l -(JQ,jh2) 

T{h) BoC(h) X b ^Y b -L{h) 



T{f 2 ) Bo£(/ 2 ) XWtoCfo) 

c{h) _ Ih (c(hi) - m,iC{h 2 )) 

A/2) n j (^(/ii)-M2,^(M) 
= ^(^i)/^ 2 )) =m( } 
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□ 

This gives the following algorithm: 

Convex bivariate decomposition 

Input: / = /i// 2 6 K(X,Y), where N(fi) C N, N(f 2 ) C N and N is normalized. 
Output: A decomposition of / if it exists, with / = uoh, u = U\/u2 with deg u > 2, 
and h = /11//12 non-composite. 



(1) Compute F = T(/i)/T(/ 2 ) 

(2) Decompose F = u(H). 

(3) Return / = where ft = ^TT^T = — ^ri'^-iYrri e K ( X > r )- 



£-\H 2 ) C£ _ 1 ( J ff 2 ).£- 1 ( J ff 2 ) 



Proposition 21. T/ie algorithm Convex bivariate decomposition is correct. 

Proof. This follows from Proposition [20] □ 

Proposition 22. The algorithm Convex bivariate decomposition uses one decompo- 
sition of a rational function of degree at most 9\N\ and O(o 2 5) bit operations. 

Proof. We apply jBLlOi Theorem 2] to N. □ 

The proof of Theorem [19] comes from Proposition [21] and Proposition [22] and 
complexity results given in Theorem [1] and |ChelO| Theorem 2]. 
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